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NOTE  ON  ASYMPTOTIC  EXPANSIONS 


Jurgen  Moser 

§1.  We  consider  a  system  or  linear  differential  equations  which  depends 
on  a  parameter  «  . 

(1)  =  Aft,  «)y  =  *  ''  {Aoft)  +  t  Aj  (t)  +  . .  }y 

wtren  y  =  (y,  , . yn  is  an  n- vector  and  A^  n  x  n  matrices.  If 

the  parameter  <  Is  considered  small  It  is  well  known  that  the  solution  of 
such  a  system  can  be  described  by  asymptotic  -  in  general  divergent  - 
seiles  expansions.  The  expansions  contain  except  for  powers  of  «  with 
rational  exponents  exponential  terms  of  the  form 

1 

a 

whore  x(t»  * )  Is  a  polynomial  In  <  ,  q  being  a  positive  Integer.  This 

theory  has  a  vast  literature  for  which  we  refer  to  the  paper  of  H.  L.  Turritln 
[lj.  While  In  many  expositions  of  tills  subject  strong  assumptions  on  the 
eigenvalues  of  Aoft)  (o.  g.  to  bo  distinct  and  different  from  zero)  limit 
the  generality  considerably,  it  is  the  aim  of  Turritln' s  paper  to  give  a 

liujUtj  t ii  *..»*■  ytiuouiJii/  vo/voiLicjiiiy  luimiiy  puuita;.  nun  icuua,  un  intz 

other  hand,  to  many  Involved  matrix  calculations  and  the  reduction  Is 
rather  complicated.  Tho  main  difficulty  seems  to  stern  from  nllpotent 
matrices  Aoft). 

In  this  note  we  study  a  question  which  is  related  to  the  presence 
of  nllpotent  matrices  Aoft)  and  which  alms  at  the  description  of  the  prin¬ 
cipal  part  of  the -exponential  8  x(t,  «)(!.-».  the  terms  of  highest  order  In 
-  iff' 

«  .  It  Is  known  that  In  case  the  eigenvalues  of  Aoft)  are  distinct  and 
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different  from  zero  that 

civ 

dT  =  Vl>*'P  +  ••• 

where  the  X  are  the  eigenvalues  of  Ao(t)  .  It  is  surprising  that  *.  , 

V  v 

which  ere  invariants  under  similarity  transformations,  play  a  role  in  the 
asymptotic  expansion  and,  thorefore,  shoulu  be  invariant  under  coordinate 
transformations 


(2) 


y  =  T{t,  t)z  . 


Such  a  transformation  maps  (1)  into 

df  = 

where 


(3)  C  -  T~  *AT  -  T  *T  . 

This  equation  (3)  defines  an  equivalence  between  A  and  B  which 

differs  from  a  similarity  transformation.  It  is  our  aim  to  characterize 

the  principal  part  of  x(l>  ')  lnvarlantl^  under  (2)  . 

q 

In  (2)  wn  admit  power  series  »  (q  -  1  ,  Z  ,  .  .  )  ao  that 

B=s  {Bo  +  •  •  } 

Is  a  series  in  fractional  powers  of  i  .  The  rational  number  m  will  be 
called  the  order  of  B  .  It  Is  our  aim  to  minimize  the  order  m  under 
transformations  (2)  and  we  define 

p-  ^ln  m  (if  positive,  otherwise  p  =  0) 

D  A 

as  the  minlrnai  values  of  m  .  That  this  minimum  Is  attained  will  be 
shown  by  a  construction  which  allows  its  computation. 

A  matrix  B  Is  called  minimal  if  m  -  p  .  It  will  be  shown  that  for 
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jj.  >  0  B  la  minimal  if,  and  only  If 

f(X,  x,  «}  --  det  (XI  -  Ba)  t  X  n 

and  that  this  polynomial  f(X,  x,  t)  In  Invariant  under  tho  equivalence  (3) 
This  result  can  bo  interpreted  In  tho  following  way:  The  equivalence  (3) 
deviates  from  a  similarity  and,  In  fact,  tho  additional  term  T  might 
Introduce  terms  of  higher  order  than  were  present  in  A  .  If,  howovcr, 

A  In  minimal,  then  theso  terma  aie  necessarily  nllpotent.  Conversely,  If 
the  highest  order  terms  in  A  are  nllpotont,  then  they  can  be  removed  and 
p  diminished. 


Wo  discuss  a  trivial  exampin:  The  matrix 
A(t,  « )  =  «' 


-P  A  -  1 

v  -  \ 


Is  nllpotont,  f,e.  all  eigenvalues  zero.  This  does  not  Imply  ihat  the 
asymptotic  expansion  la  free  from  exponential  terms  {l,e.  x  3  0)  and,  In 


fact,  one  finds  a  so'utfon 

-ft 


exp  (±«  )  ,  y2  -  (t« 

E 


1 


-S> 

1)  exp  (±«  ) 


This  shows  that  p  =  2  which  can  be  easily  obtained  by  the  method  to  be 


discussed. 

The  following  considerations  are  algebraic  In  nature  and,  therefore, 
all  series  can  be  considered  convergent.  This  result  can  be  used  in  a 
general  theory  of  asymptotic  expansions  which  will  not  be  done  here,  since 
a  treatment  of  this  kind  can  bo  found  in  Hukuhara's  paper  [2]*  . 


*  1  am  grateful  to  Professor  Slbuya  for  this  reference. 
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§2  Necessary  and  sufficient  conditions  for  minimal  B  ■ 

To  simplify  the  following  considerations,  we  assume  all  elements 
of  the  matrices  A^(t)  to  be  ineromorpMc  functions  of  t  *  .  The  complex 
variable  ranges  over  a  fixed  domain  D  . 

q 

We  denote  by  Fo  the  space  of  functions 

v 

00  — 

f(t,«)  --  Q(t)'1  'l  Mt)<q 

v-u 

whore  ,  Q{t)  are  analytic  in  T>  and  the  series  converges  in  , 

t  (-  D  .  These  functions  form  a  ring  and  quotients  of  such  functions  will 
have  'he  form 

v 

oo  — 

-  1  \  ■  q 

f  :  Q  l  Mt)< 

l/  =  “V0 

which  form;.  1  he  space  Tq  .  Siniilaily,  the  matrices  T  ,  whose  elerr.er.ts 
!iit  in  Fo  q  ,  1^  ,  define  the  spaces  rospectl  J'-\y .  The- union 

of  ail  M'1  for  a.  I  ,  i  ,  .  .  will  bo  denoted  by  M  . 

The  Important  oroperty  of  f  1;;  that  addition,  subtraction,  division 
and  differentiation  doe:;  not  lead  out  of  l,t(  . 

A  matrix  T  (-.  My  .  i.e. 


00 

..  Y  nvT 
-  Lj  n  Tv 
v-O 


is  cfllod  a  unit  if  T  *  also  lies  in  W  ^ 

o 

If, 


This  is  the  case  if,  and  only 


dot  To  (t)  0  In  D  . 


*  Under  this  assumption  It  Is  not  necessary  to  shrink  the  t-domaln  which 
is  usually  done  in  order  to  avoid  that  eigenvalues  cross  each  other. 
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Lemma  1:  If  det  T  ^  0  Is  an  element  of  M  ,  then  It  can  oe  written  In 
the  form 

(4)  T  =  Pit,  i)  caQ(t,«) 

where  P  ,  Q  are  units  In  (l.e.  det  Pa  ♦  0  ,  detQsipO}  and 

«  a  =  dlag  («  01  ,  , ai  ,  . , .  ,  t an)  (aj  <  aj  <  . . .  <  a  ) 

and  qa  are  integers. 
v 

Proof;  Without  loss  of  generality  -  can  assume  q  =  1  ;  otherwise 

\ 

q 

let  rj  a  i 

The  matrix  Q  will  be  built  up  from  “elementary  matrices"  consist¬ 
ing  of  permutation  matrices  (independent  of  t,  <)  whOBe  determinant  is 
*1  and  matrices  which  differ  from  the  unit  matrix  In  one  off  diagonal  ele¬ 
ment,  All  these  matrices  us  well  as  their  product  are  units,  since  their 
determinant  is  *1  , 

Nov/  determine  at  as  the  greatest  Integer,  Buch  that  <  a'  T  is  a 
power  series  In  <  ,  1.3.  i  a,T  t  M’c  .  Thus  the  constant  term  in 

—  a  i 

t  T  is  non-zero.  Applying  a  permutation  matrix  Q  ,  we  can  achieve 


that 


TQj  —  <  (s,  ,  sj  ,  ...  ,  sn) 

where  =  s^(t,  *)  are  c°lumn  vectors  containing  no  negative  powers 
of  i  and  s,  (t,  0)^0  . 

Now  choose  a*  as  the  maximal  Integer  such  that 


TO  =  (‘ Q* r  ,  «  a* r.-,  , 


V 


where  Q  Is  any  product  of  elementary  matrices  and  c  are  column  vec- 

V 

tors.  Obviously,  at  >  .  After  having  defined  ij  ,  . . .  ,  we 
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maximize  under  all  products  of  elementary  matrices  Q  such  that 

a.  i  a«  a,  Q. 

,,,  ,  Oj  1-1  1  I  I  , 

(5)  Tw  =  («  ‘p,  .  ,  «  P|_j  *  Pi  ,  *  P1+]  ,  ...  ,  t  pn) 

and  the  elements  of  .he so  column  vectors  belong  to  F0  *  . 

The  existence  of  such  a  maximal  n^  follows  from  the  assumption 

det  T  =  c(t)«Y  +  ...  c(t)  *  0  . 

Y  gives  an  upper  bound  for  since 

«1  +  *i  +•••■*'  *1  ^  +  (n  -  l  +  1 )  <  y  . 

One  proves  by  induction  that  the  rank  of  the  matrix 

(Pi(t,0)  ,  ....  ,  pn(t,0)) 

over  the  mcromorphlc  functions  is  at  least  t  .  Otherwise  one  can 

achieve  by  an  elementary  matrix  that  p  0)  s  o  for  some  v*  <  I 

which  means  that  can  be  enlarged.  Thus  foi  =  n  one  has 

TO  =  (p,  ,  - ,  P  )  ta  -  P(t,  «)  tn 

n 

det  P(t,  0)  *  0  . 

Observe  that  with  Q  also  Q  1  is  an  elementary  matrix,  hence  a  unit. 
This  proves  the  Lemma. 

Remark:  The  above  representation  is  not  unique  but  It  Is  easily  seen  that 

the  a  are  unique. 
v 

Theorem  1:  Let 

A  -  t'P  {Ao  (t)  +  ....  } 

{6>  1 

^  B  =  «  P  {Bo<0  +  ■  •  •  •  } 

be  two  matrices  in  M  which  are  equivalent,  i.e.  related  by  (3)  with  a 
T(E  M  ,  p  Is  assumed  to  be  a  positive  rational  number. 
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det  (XI  -  Ao  (t)}  =  del  (XI  -  B0  (t)) 

q  i  q« 

Proof:  Assume  A  contains  integer  powers  of  «  and  B  in  <  then 


define  tj  by 


«  =  n 


so  that  A  ,  B  contain  only  Integer  powers  of  q  .  Let  T  be  the  trans¬ 
formation  relating  A  and  B  by  (3)  .  Represent  T  in  the  form  (5) 


T  »  Pt|  <5 


and  lot 


A  =  P*  AP*  P’  P 
8  =  Q  B  O"*  +  Oq"* 

then  &  is  equivalent  to  A  and  B  equivalent  to  B  .  Moreover,  since 
P  ,  Q  are  units,  P  *  ,  Q  are  power  series  and  so  arc  the  terms 

P  *P  and  Q  Q  *  •  Since  p  »'  o  it  follows  that 
A  ~  f  ^  {  Pp  Ao  Pq  r  .  .  .  } 

B  -  «  t  Qo  Bo  Oc>  +  -  •  • } 


/dot  (XI  -  Ao)  =  dot  (XI  -  Ao) 

(?)  < 

^  dot  (XI  -  B0)  =  det  (XI  -  B*  1  . 

Finally,  A  ,  B  are  related  by  the  similarity  transformation 

o  -n  t  a 
B  =  rj  n  t] 

from  which  it  follows  that  the  first  determinant  in  (7) 


det  (XI  -  A)  =  det  (XI  -  Ac) 
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agrees  with  the  second,  which  proves  the  Theorem  1  . 
Definition:  A  matrix 

A  =  i'p  {Ao  +  . ..  }£  M  ;  p  >  0 
is  called  minimal.  If  for  each  T  £  M 

B  =  T_1AT-T=1T  =  {B0+..} 

one  has 


m  2.  p  >  0  . 

Otherwise  B  Is  called  not  minimal,  for  p  >  0  .  For  p  >  0  A  Is  defined  to 
be  minimal  unconditionally. 

Consoquence  of  Theorem  Iz  11  A  is  not  minimal,  then  Ao  Is  nilpotent. 
Proof:  Since  A  Is  not  minimal,  there  is  a  T  <E  M  such  that 
B  =  T_  1  A  T  -  T~ 1  T 

does  not  contain  terms  of  order  >  m  in  «  1  ,  1.  e.  in  (6)  one  has 

Bo  a  0  .  Hence  by  the  Theorem  1 

det  (M  -  Ao)  =  det  (X.I  -  0)  =  X.n  , 
which  expresses  that  Ao  Is  nilpotent. 

Theorem  2;  Ever/  matrix 

A  =  ."P  {Ao (t)  ...}  £  M 
Is  equivalent  to  a  minimal  matrix  B  . 

Proof;  This  theorem  requires  a  construction  and  guarantees  the  existence 
of  a  minimal  value  p  of  m  .  We  make  use  of  a  procedure  which  has 
been  used  by  Hukuhara  (2]  in  a  similar  connection.  Ths  use  of  Theorem  ! 
will  simplify  Hukuhara's  Lemma. 


110.1 
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It  is  well  known  that  each  component  of  y  satisfies  a  single  dif¬ 
ferential  equation  of  order  <_n  ,  to  construct  such  a  differential  equation 
for,  say,  y,  ,  we  set 


Zi  =  Yi 


*»  =  dh  =  £  av  <‘.‘)yv 


V=1 

and  expross  yj  ,  y2  In  terms  cf  z-t  ,  z2  if  aJ2  t  U  .  If  a;  *  ■*  0  for 
some  v*>  1  lenumber  y2  ,  ,  y  ,  to  reduce  this  to  the  previous  case. 

If  a^  b  0  lor  v  -  2  n,  Zj  =  =  a,  ,  y,  which  represents 


the  desired  differential  equation.  Otherwise  let 

dz.  \  ■  - 
*  dt  - 

which  is  a  linear  combination  of  y. 


y  +  y 

V  V  LJ  tv  X, 


Yl 


combination  of  zt  ,  z{  ,  y3  ,  ,  y 


y  and  hence  a  linear 
n 


If  the  coefficients  of  y3  , 


.  .  .  ,  yn  are  not  all  zero,  wo  can  assume  by  renumbering  that  the  coefficient 
of  y5  Is  not  Identically  zero,  which  allows  to  express  y3  ,  y2  ,  y3  in 


terms  of  z3  ,  z2  ,  Zj  ,  z,  ,  .  .  ,  yR  . 


dz. 


This  proceos  of  elimination  can  be  continued  until  — r-T  is  a  linear 

dt 


combination  of  Z|  ,  .  .  , 
z  =  Bz 

whei  e 


.  If  t  =  n  ,  one  obtains  for  z  a  system 


(8) 


B  = 


1 

b, 
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If,  hovrever,  r  <  n  one  can  continue  the  same  procedure  by  defining 


T  +  l 


T  +  l 


introducing 


dz 


T  4-1 


and  eliminating  the  further  y^  for  v  >  r  .  This  leads  to  a  system 
i  =  8{t,  < )  z 


whore 


(9) 


B  - 


Bn  0 


Kl 

Herf!  the  8^  (v  l  ,  k)  aro  Square  matrices  of  t  ho  type  (  H) 
While  tlie  (v  --  i»)  contain  nonzero  oicmont;;  only  in  the  last  column. 

The  variables  y  unci  z  are  related  by  linear  equations,  the 
coefficients  of  winch  are  obtained  from  those  A  by  addition,  multiplication, 
division  and  differentiation,  hence  belong  to  f  “  li  A  (y  m"  .  Since  the 
relation  is  invertible,  we  have 
y  =  Tz 


with 

0°)  T  (-  Mq  ;  del  T  *  0 

1.  e.  A  and  B  are  equivalent. 

It  Is  easily  shown  that  B  can  be  transformed  Into  a  minimal  matrix 
by  a  diagonal  matrix  of  the  form 

E  =  <‘Y*  ,  ...  ,  <Vn  )  . 


#104 
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Flrst.  one  can  easily  achieve  that  B  for  v  >  p  does  not  contain 
negative  powers  of  «  by  a  transformation 
<1 ' )  E  =  dlag  (E,,  ,  E« . E^J 

where  the  E  are  square  matrices  (corresponding  to  the  block  represents- 

V  V 


tlon  (9)  of  B)  and 

E 

V  V 


-p 

V 


-  1  ry  'V 

The  matrix  E  B  E  Is  obtained  from  B  by  replacing  B  by  B  i  . 

vp  vp 

Choosing  for  |t  an  appropriate  increasing  sequence,  one  can  achieve  that 


B  (-  Mo  • 
vp 

If  In  (9)  no  8^  contains  negative  powers  of  c  ,  then  B  (-  M^ 
and  B  is  minimal  by  definition.  Otherwise,  If  B  ,  which  Is  of  the 

vv 

form  (8)  ,  cor  tains  negative  powers  of  «  ,  determine  the  smallest 
rational  number  y  such  that 


bj 


bi 


z 

Yv 

c 


b_  t 


Vv 


contain  no  negative  powers  of  t  .  Thus  at  least  one  of  these  terms  does 


not  vanish  for  r  =  0  ,  Fly  definition,  v  s  0  is  a  rational  number.  Let 

*  '  V 


y  =  Max  y 

V  V 


and  apply  the  transformation  with 


E 


-Y 

dlag  (1  ,  t  , 


Since 


V  V 
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0  1 


vv  vv 


with 

c  =  b 

and  at  least  one  of  them  does  not  vanish  at  <  =  0  .  Since  the  character¬ 
istic  polynomial  of  .  E  B  E  la  XT-Cj,  XT  1  c  It  Is 

VV  VV  VV  r 

not  Identically  \T  for  <  =  0  and  for  some  v  .  Thus  we  obtain 
E*  1  B  E  =  «~V  {Co  + . } 

when 

dot  (XI  -Cp)  $_  \n 

which  proves  that  E  1  B  £  is  minimal.  From  theorem  1,  It  follows  that 
Y  >  p  .  This  proves  theorem  2  . 

Definition;  For  a  given  mat-lx  A  let 

B  =  «"m  {B0  + - } 


be  an  equivalent  minimal  one.  Then  define 

...  |  m  if  m  >  0 

^  =  |1|A)  =  |  0  If  m<0  . 


The  number  p  has  by  definition  an  Invariant  meaning.  In  fact,  it 


is  the  minimal  order  of  matrices  which  are  equivalent  to  A  . 
quence  of  the  proof  of  theorem  2,  we  show 


As  a  conse- 
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Theorem  3:  For  a  given  (not  necessarily  minimal)  matrix  A  fc  .»*  one  can 
compute  (r  =  p(A)  by  transforming  it  into  the  form  (9)  by  a  transforma¬ 
tion  T  £  and  determining  the  smallest  m  >  0  such  that 

m  _ 

riet  rxi  -  •  n\C  xa. 

Then  m  =  h-  • 

Proof:  Since 

C  =  1  B  E  =  i_|1  (Gi  +  . . .) 

is  minimal  It  follows  that 

x"  «  det  (XI  -  Co )  =  det  (XI  -  il)  |  (  __  Q 

which  implies  m  =  p.  . 

Theorem  1  Implies  that  the  characteristic  polynomial  of  a  minimal 
matrix  Is  Independent  of  the  representation  If  p  >  0  ,  which  proves  the 
invariance  of  the  nonzero  eigenvalues  of  a  minimal  matrix  under  equiva¬ 
lence. 


dt 


It  is  easily  verified  that  these  nonzero  eigenvalues  agree  with 
where  x  was  defined  in  the  introduction.  If  p  =  0 


«  =  0 

the  eigenvalues  do  not  play  any  role,  in  fact,  in  that  case  A  is  equiva¬ 
lent  to  B  a  o  .  To  prove  this  statement  let 

il L_ 


dt 


=  Ay 


where  A  £  Ms  does  not  contain  negative  powers  of  t  .  It  Is  known 


that  In  this  case  the  matrix  solution  of 


4  Y  =  A  Y  ;  Y  =  I  for 
at 


0 

=  1  . 


Ji i 


belongs  to  Mg 


and  is  a  unit,  since  det  Yo 
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If  one  makes  the  transformation 


T  =  Y 

y  =  Tz  =  Yz 

every  solution  y  goes  over  into  a  vector  z  which  Is  independent  of  t  > 
namely 

Ay  =  y  -  Yz  f  Yz  =  AYz  +  Yz 


hence 


z  =  0 

which  proves  the  statement 

§3  We  mention  without  proof  that  one  can  obiain  the  principal  part  of 
those  exponentials  y  which  start  with  terms  of  lower  order  than  s  11 
by  refining  the  above  method.  Instead  of  minimizing  the  order  of  A 
as  defined  above,  one  can  minimize  a  “matrix  order"  or  order  which  is 
defined  as  follows:  Assuming  for  simplicity  det  A  *  0  wc  can  represent 
A  by  the  previous  Lemma  as 


where  m  -•  diag  (mj  ,  m*  ,  .  =  ,  m  )  m,  >  mj  .  Ordering 

m  lexicographically,  i.e. 

m  ;>  ifi  ifm  =  rn  u  <  k 

V  w 

and  m  >  m 

K  K 

one  can  find  conditions  on  minimal  matrix  m  which  lead  to  a  description 
the  principal  parts  of  the  other  exponentials. 
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